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Abstract : 

We use improved lattice actions for glue, light quarks and heavy quarks for which 
we use lattice NRQCD to compute hadron masses. Our results are in good agreement 
with experiment, except for charmed hadrons. It seems that charmed quarks are not well 
approximated as heavy quarks nor as light quarks. 



PACS numbers: ll.lSHa, 12.38G, 14.40G and 14.40Gx 



1. Introduction 

A notorious feature of the perturbative theory of lattice QCD is that the perturbation 
series converges slowly at a moderate lattice spacing. Predictions based only on first two 
or three terms in the series are in most cases still not reliable. A recent study[l] shows that 
a partial reason for this is the large contribution from tadpole diagrams, which are absent 
if one uses dimensional regularization instead of lattice regularization. Based on this the 
method of tadpole improvement [1] for lattice perturbation series has been suggested. This 
has many practical implications for studies of lattice QCD. Most importantly it makes the 
idea of improved actions practicable in simulations of lattice QCD. The idea of improved 
actions was proposed a long time ago [2], to systematically improve lattice actions using 
perturbation theory to remove the effect introduced by the finite lattice spacing a. If such 
an improvement really works, one can simulate lattice QCD on a coarse lattice with a 
lower cost in CPU time, but still obtain reliable results for physics in the continuum limit. 
Several simulations [3-6], where the tadpole improvement is implemented, already show 
that this is possible. 

With an improved action for glue at the one-loop level and the tree-level improved 
action for NRQCD the mass spectrum of charmonium is obtained in [3] from lattices with 
a lattice spacing a between 0.5fm and 0.2fm, and it agrees with experiment. Similar results 
were also obtained in other simulations [4], where the mass spectra of bottonium and the 
Be meson, in addition to that of charmonium, were obtained. With an improved action 
for light quarks, the D234 action, the mass spectrum of light hadrons is calculated from 
simulations on coarse lattices [5], and good agreement with experiment is found. In general 
improved actions are not unique, in [6] another improved action for light fermion, which 
was proposed long time ago [7], was employed, the mass spectrum of light hadrons was also 
successfully calculated[6] , where the gluonic action improved at tree- level was used. Unlike 
the widely used SW-action for light quarks [8], in which the effect of 0(a) is removed at 
tree- level, the actions used in these simulations are improved up to O(a^) at tree- level. 

In this work we will present our results for the mass spectrum of light hadrons, strange 
hadrons, and hadrons containing a b- or c-quark. In calculating the mass spectrum we use 
the improved action for glue at the one-loop level and the action of [7] for light quarks. 
For the heavy b- and c-quarks we employ the action of lattice NRQCD. The actions for 
quarks are all improved at tree-level upto O(a^). We use the quenched approximation 
for light quarks. Two lattices are used with the size 6^ x 12 and 8^ x 16. Their lattice 
spacings determined by the charmonium spectrum are 0.41fm and 0.23fm respectively. 
The whole calculation is performed on two UNIX workstations, which have computational 
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power equivalent to that of a high-performance personal computer (like a Pentium Pro 
200). Our paper is organized as follows: In Sect. 2 we introduce the actions used in this 
work. Our results for light hadrons and strange hadrons are given in Sect. 3. The results 
of c-flavored and b-flavored hadrons are given in Sect. 4. Sect. 5 is the summary. 



2. The Improved Actions 



2.1 The gluonic action: 

We take the one-loop improved action for glue [9], where the action consists of pla- 
quette, rectangle and paralellogram terms and is accurate up to errors of 0{ala^,a'^). 
Implementing tadpole improvement the action becomes [3] 

S{U) = (3Y1 ^R-eIV(l - Upi) + Prt ^R-eIV(l - Utr) + Ppg ^R-eIV(l - Upg), 

pi rt pg 

(^rt = -7^(1 + 0.4805a,), (5pg = -4o.03325a„ (2-1) 
1 , ,,1 ln(iReTr(t/pi)) 

We used this action to generate gluonic configurations at /? = 7.4, on a lattice whose 
size is 8^ x 16, and also at /? = 6.8 with the size 6^ x 12. The parameter uq is determined 
by self-consistency. It is 0.8631 and 0.8267 at = 7.4 and /? = 6.8 respectively. The 
pseudo heat bath method[10] was used to update the links, and the three SU (2) subgroups 
were updated 3 times in each overall update step. For each lattice we generated 100 
configurations for our mass calculations. 

2.2 The improved action for light quarks: 

The improved action proposed in [7] can be written as: 



-Slight = - E I ™V'(^)V'(^) + - ci/^Tm^) + r Y 

(2.2) 



(2) 

where and Ajj, are lattice derivatives with the gauge link Un{x): 

^mV'(^) = \iU^{x)ij{x + fi)- Ulix - fi)i;{x - fi)), 

Al^^{x) = U^{x)'4){x + A) + Ul{x - il)iij{x -fi)- 2ii;{x). 

The parameter ci is determined to be 1/6 at tree-level to remove the effect at O(a^), m 
is the mass parameter for a quark. The last term in the action with the parameter r 
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is introduced in analogy to the Wilson term in the Wilson action to solve the doubling 
problem of lattice fermions. However, the doubling problem is not totally solved. In the 
free case one can solve the equation of motion determined by the action to see whether the 
doublers are removed or not. An analysis in [8], and the analysis in [5] for the D234 action, 
show that in the low-energy regime the actions describe one particle in the sense that the 
propagator has only one pole. In high-energy regime there are additional "unphysical" 
poles. As we are only interested in the low-energy regime as we are in this work, we 
can expect that the effect from other "unphysical poles" to be negligible. We will take 
r = 1/6. With this choice the action is the same as employed in [6] , in which it is shown that 
there is indeed no effect in the low-energy regime which can be related to the "unphysical 
poles". In this work we used the stabilized biconjugate gradient algorithm[ll] to calculate 
propagators for light quarks. On our lattices this algorithm is at least three times faster 
than the conventional conjugate gradient algorithm. 

2.3 The Action of Lattice NRQCD 

Heavy quarks whose mass is larger than 1 in lattice units cannot be simulated directly 
as above with reliable results. To simulate them one uses the heavy quark effective theory, 
HQET. The formulation of HQET for hadrons with zero velocity is equivalent to that 
of NRQCD on the lattice, except that the expansion parameters are different. As we 
will only create hadronic states on lattice with zero space-momenta, we may use lattice 
NRQCD for heavy quarks like the b- and c-quarks. On the lattice one needs to calculate 
the propagator of heavy quarks satisfying its appropriate evolution equation. We take the 
evolution equation proposed in [12,13]. The propagator G{t) (where G{t) = for t < 0) 
can be calculated on the lattice as: 



TTif, 



n 



where 



G{t + l)^{l-^) Ul(l-^] {l-6H)G{t), 



Ha 



2n J ^ V 2n 

A(2) 



(2.4) 



2Mq 



2Mq 24Mq 16nM2 



(2.5) 



In (2.5) A^^) is the lattice version of the continuum operator Df. Mq is the mass 
parameter for the heavy quark Q. The last two terms in 5H are the correction terms to 
remove the effect at order O(a^). The first term is responsible for spin-splitting in the 
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mass spectrum, where B is the chorommagnetic field. We use definition of B in terms of 
gauge hnks given in [12], and it is also improved up to errors of order 0(a^,^^a^). The 
parameter n is introduced to avoid numerical instability when high-momentum modes 
occur. With propagators calculated with Eq.(2.4) we reach an accuracy of order 
in the mass spectrum. With lattice NRQCD we have calculated the mass spectrum of 
quarkonium at the above /?-values[4], and have determined that at /? = 7.4 the mass 
parameters for b- and c-quarks are = 4.6 and Mc = 1.4 respectively. We will use these 
parameters for our calculations of mas spectra of b- and c-flavored hadrons. We take n = 3 
for c-quark and n = 1 for b-quark. 

In our calculations of quark propagators with the action in Eq.(2.2) or with Eq.(2.4) 
tadpole improvement is implemented. 

3 The Mass of Light Hadrons 

We create hadronic states on the lattice by using standard local operators Oh{x): 

0>^{x) =u{x)'y''d{x), (3.1) 
Op{x) =eabcUa{x) [uj {x)C'y5dc{x)] 

for TT, p and proton respectively. In Eq.(3.1) u{x) and d{x) stand for u- and d-quarks, the 
indices a, b and c are color indices. With these operators one can measure the corresponding 
hadron correlations Cnit), in which the hadron has zero space- momentum and its spin is 
averaged if it has spin. We only use local sources to calculate Cnit), and u- and d-quarks 
are taken to be degenerate in mass. The hadron correlation Cnit) is fitted in a certain 
time-interval as: 

an (e-"""* + e'^^^^^-^A for H = tt, p, 

T (3-2) 
a^e with t < for H = proton, 

where T is the lattice size in the time-direction. 

Using the lattice configurations on an 8^ x 16 lattice for P = 7.4 we calculated these 
hadron correlations varying the mass parameter for light quark from —0.5 to —0.75, and 
fit them according to Eq.(3.2). The negative values of the mass parameters simply reflect 
the difference in sign between the mass parameter and the "Wilson" parameter r. The 
range of light quark mass parameters investigated here are the same order as the mass of 
the s-quark. The fitting window for C.j^{t) and Cp{t) is chosen from t = 6 to t = 10, for 
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the proton we can see a plateau in the region of t = 4 to t = 8, so we take the window to 
be from t = 4 to t = 7. As examples which show our fit, we plot the measured Cnit) for tt 
and p with the fitted results in Fig.lA and Fig. IB respectively. Fig.lC is the plot for the 
effective mass of the proton. In all these figures the results are at m = —0.65. Our results 
of fitted masses are given in Table 1. 



Table 1 



/? = 7.4 


m= -0.5 


m = —0.6 


m = —0.65 


m= -0.7 


m = —0.75 




1.142(9) 


0.941(11) 


0.832(14) 


0.713(18) 


0.577(27) 


rup 


1.362(20) 


1.211(24) 


1.138(32) 


1.068(46) 


1.010(95) 


mp 


2.23(9) 


1.96(7) 


1.82(6) 


1.69(7) 


1.56(16) 



We also constructed hadron correlations for pion with a minimal lattice momentum 
|p| — T"' where L is the lattice size in the space-direction. Measuring these correlations 
one can investigate the dispersion relation -E^(p) — = c^|pp. We find that c is close 
to 1. For example, we obtained c = 1.02(7) at m = —0.5 and c = 1.00(18) at m = —0.75. 
This fact indicates that rotational invariance is restored on the lattice. To find the masses 
in the chiral limit, we fit the hadron masses with: 

mp = vo + viml, mp = uq + uiml + U2ml. (3.3) 

We find that the masses given in Table 1 are well described by these relations as shown 
in Fig.2. For the mass of the proton we clearly see the eff'ect from the cubic term. With 
these fits one can extrapolate the masses to the chiral limit. 

If we assume that the quark mass is obtained through multiplicative and additive 
renormalizations of the mass parameter, then the square of m^r should be linear in the 
mass parameter m. With the data in Table 1 we find good agreement with this relation. 
Hence we are able to determine the critical mass-parameter mo at which itIt^ = 0. We 
also determine the mass parameter of the strange quark from the experimental value 
m\-/inn?p = 0.412. These parameters at /? = 7.4 are: 

mo = -0.830(40), m^ = -0.6658. (3.4) 

With the nis above we calculated the propagator of the s-quark and then the hadron 
correlations for K, K* and ^, where Cxit) and CK*{t) are calculated with these mass 
parameters in Table 1 for light quark. The masses obtained for niK, ttik* with non-zero 
mass for the u-and d-quark are given in Table 2. 
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Table 2 



(3 = 7.4 


m = -0.5 


m = -0.6 


m = -0.65 


m= -0.7 


m = -0.75 


rriK 


0.979(12) 


0.870(13) 


0.814(14) 


0.755(16) 


0.694(20) 


rriK* 


1.239(25) 


1.163(29) 


1.127(33) 


1.092(40) 


1.061(50) 



We fit these masses with: 

= bo + bi{m — mo), uik* = cq + ci{m — mo) (3-5) 

to extrapolate to the hmit of zero-mass u- and (i-quarks. The data and the fitted hne 
are drawn in Fig. 3. The data points are well described by the linear fitted line giving 
confidence in the extrapolation. 

Doing the same calculations of Cnit) with our configurations at /3 = 6.8, for which 
we choose the region of the light quark mass parameter m from —0.8 to —1.1, we obtain 
hadron masses which are given in Table 3. 



Table 3. 



P - 6.8 


m = -0.8 


m = -0.9 


m = —1.0 


m = —1.05 


m = —1.1 




1.353(9) 


1.179(10) 


0.986(11) 


0.877(13) 


0.754(14) 


nip 


1.703(21) 


1.587(29) 


1.470(48) 


1.409(62) 


1.340(81) 


rriN 


2.770(77) 


2.523(77) 


2.262(79) 


2.113(86) 


1.941(106) 


rriK 


1.218(10) 


1.125(10) 


1.028(11) 


0.977(12) 


0.924(12) 


rriK* 


1.614(27) 


1.554(33) 


1.495(43) 


1.467(49) 


1.439(56) 



At this /?-value we have determined: 

mo = -1.233(56), m^ = -0.9513. (3.6) 

With our results we obtain hadron masses in the chiral limit and the mass m,^ which 
depends only on mg. They are given in Table 4 as a ratio to m^, which is itself given in 
lattice units. 



Table 4 





P = 6.8 


/? = 7.4 


Exp. 


arup 


1.198(67) 


0.885(45) 




mp/nip 


1.16(10) 


1.45(11) 


1.22 


rriK/mp 


0.63(4) 


0.64(3) 


0.64 


rriK* /TTip 


1.14(7) 


1.13(6) 


1.16 


m^/rup 


1.27(8) 


1.26(9) 


1.32 
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Having chosen the parameter to fit the mass of the our predictions for the meson 
masses are in good agreement with experiment, although the quenched approximation was 
used. The mass of proton obtained at /5 = 7.4 is 20% larger than the experimental value, 
while that from j3 = 6.8 agrees with the experimental value within the statistical error. 
A possible reason for the deviation at /5 = 7.4, apart from the quenched approximation, 
could be the effect of the finite volume. Another possible reason is the use of a local source 
in our work, a result of which is that the signal of the ground state in Cp{t) is rather weak, 
and the fitting quality of Cp{t) is worse than that of meson propagators. This situation 
may be improved if smeared sources are used. 

From the calculated p-mass in lattice units we can determine the physical lattice 
spacing. We obtain: 

a-^ = 0.64(4)GeV, at = 6.8, 

, (3-7) 
a-i = 0.87(4)GeV, at /? = 7.4. 

It is instructive to compare the lattice spacings determined from the quarkonium system. 
In our previous work [4] we obtained: 

a-}^ = 0.477(5)GeV at /? = 6.8, 

(3 8) 

a'!^ = 0.749(4)GeV, a^^ = 0.861(5)GeV at /3 = 7.4. 

The lattice spacing at /? = 7.4 determined from p is close to that obtained from T, while 
there is a deviation of 14% in comparison with aj^^. At /? = 6.8 the deviation is at 
30%. Similar lattice spacings were also obtained with D234 action in [5]. Possible reasons 
for these deviations include: the effects neglected in the quenched approximation being 
different in different systems, the effect from higher orders in lattice NRQCD can be large 
for charmonium system because the charm mass is not enough large. In general one should 
not be surprised that the lattice spacings determined from different systems are different 
when approximations are used in the calculations. 

From the determined dependence of rrip on ttIt^ we can obtain the quantity J = 
mpdm^/dm^ at the experimental value mp/rriT^ = 1.8. This quantity is introduced in [14] 
to judge the quality of the quenched approximation. For the "real world" it is J = 0.48(2), 
while the "world averaged" result from the quenched approximation is J = 0.37. From 
our results we obtain: 

J = 0.375(87) at /? = 6.8, J = 0.363(63) at /5 = 7.4, (3.9) 

which are in consistent with J = 0.37. 
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4. The Mass Spectrum of c- and b- flavored Hadrons 



In this section we will present only our results at /3 — 7.4 in detail. At P — 6.8 there 
is a problem for simulations of 6-quark with lattice NRQCD. We find large correlation 
lengths, so that a significant finite- volume effect should be expected. On this lattice we 
calculated only masses of c-flavored hadrons and we will mention these results briefly. With 
the configurations at /? = 7.4 we will calculate masses of the hadrons S*, Bg, B* and 
Afo and masses of the D*, Dg, D* and Ag. 

For the hadrons we use operators analogous to the operators listed in Eq.(3.1) to cre- 
ate these hadronic states on lattice. The difference is that a light quark field in Eq.(3.1) is 
replaced by a heavy quark field Q{x) and the field Q{x) only has two nonzero components. 
The propagation of the field is determined by Eq.(2.4). For light quarks we use the prop- 
agators of the last section. For calculating propagators of heavy quarks we use a source 
smeared by a gaussian: 



1 + 



(2) 
i 



(x,y,t). (4.1) 



The smearing parameter Us is fixed at = 10, while the parameter e is adjusted so that 
the smearing radius is about the half of the simulated system. Using a smeared source is 
essential for our simulations. With a local source for heavy quark propagators the signal 
of hadron correlations is overwhelmed by the statistical noise after 3 or 4 time slices from 
the source. We only calculate hadron correlations with smeared sources, while the sinks 
remain unsmeared. The measured hadron propagators are fitted to the form 

Cnit) ~ ane-^"', (4.2) 

where the index H refers to the hadrons mentioned above. With the fitted energy Eh one 
can obtain the hadron mass Mh 

Mh = Aq+Eh, (4.3) 

where Ag is the difference between the renormalized mass of the heavy quark and the zero 
point energy of the heavy quark on lattice, and it does not depend on the type of hadron. 
Therefore one can predict the mass difference between two types of hadrons if the E'h's 
are known. 

Although we used a smeared source for heavy quark propagator, the signal is buried in 
statistical noise after 10 or 11 time slices from the source for mesons and after 8 or 9 time 
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slices for baryons. For the meson we already see a plateau in the effective mass. We will 
take the fit- window for meson propagators from 6 to 9. For baryons we take the fit-window 
from 5 to 7. In Fig.4a, Fig. 4b and Fig. 4c we give our effective mass plots for B, B* and 
Afe respectively. These plots are for propagators at the mass parameter m = —0.65. Our 
results for Eh of various hadrons with nonzero mass of light quarks are given in Table 5. 



Table 5 





m= -0.5 


m = —0.6 


m = —0.65 


m= -0.7 


m = —0.75 


Edo 


1.069(10) 


0.990(10) 


0.950(9) 


0.909(9) 


0.863(9) 


Ed* 


1.143(8) 


1.067(9) 


1.030(9) 


0.991(10) 


0.948(10) 


Ea. 


1.976(16) 


1.810(17) 


1.727(17) 


1.643(15) 


1.561(17) 


Eb" 


1.043(14) 


0.967(12) 


0.928(12) 


0.894(12) 


0.856(13) 


Eb 


1.066(11) 


0.999(9) 


0.962(9) 


0.927(9) 


0.891(10) 


EAb 


1.892(12) 


1.713(12) 


1.617(14) 


1.504(16) 


1.381(22) 



For c- and b-flavored hadrons with s-quarks we have: 

Eb^ = 0.920(11), Eb* = 0.952(9), 

(4.4) 

Ed, = 0.937(9), Ed* = 1.018(9). 

To extrapolate to the zero light quark mass we assume the dependence of Eh on the mass 
parameter m to be: 

Eh = Ef + E^i\m - mo), (4.5) 

and use this relation to fit our data. In Fig. 5 A and Fig.5B we show the fits to the data 
of Table 5. Except for A;,, our results in Table 5 can be well fitted with this relation. In 
Fig.5A the fitted line for A;, is determined only by the first 4 data points. The fit has 
= 0.76. If we include the last point into the fit, the jumps to 3.39, and E^^^ changes 
from 1.218(7) to 1.250(9). A reason for this may be the effect of higher orders of (m — mo) 
omitted in Eq.(4.5). In our final results we will use the fit with the 4 data point for A^. 
We obtain the energies at zero mass of light quarks: 

4°) = 0.796(4), 4°2 = 0.836(3), 4°^ = 1.218(7), ^^^^ 
E^^l = 0.800(3), 4* = 0.889(3), 4°, = 1-428(5). 

It is interesting to note that the energies given in Eq.(4.4) and in Eq.(4.6) for b-flavored 
hadrons are close to those of c-flavored hadrons. These energies Eh^ can be expanded in 
the inverse of Mq with a leading order Mq and are calculated with an accuracy of Mq^ 
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in this work. The above fact indicates that the effect from the next-to-leading order and 
from higher orders is small in these spin averaged masses. However, as we will see, this 
is not true for spin-splittings. With the results in Eq.(4.4) and in Eq.(4.6) we are able to 
predict the mass differences. For doing this we take the lattice spacing determined by the 
bottonium system for b-flavored hadrons, 

a^^ = 0.86GeV. (4.7) 

This spacing is the same within errors as the lattice spacing determined from nip in the 
last section. We obtain 

Mb* - Mbo = 34(6)MeV, Mb - M^o = 107(13)MeV, 

(4.8) 

Mb* - Mb, = 27(17)MeV, Ma, - M^o = 363(9)MeV. 
These results should be compared with the experimental results 

Mb* - Mbo = 46MeV, Mb - Mbo = 91MeV, 

(4.9) 

Mb* - Mb, = 47MeV, Ma, - Mgo = 363MeV. 

We find that our results for Ma, — M^o and Mb^ — Mbo agree well with experimental 
results, while the spin-splitting Mb* —Mbo and Mb* —Mb, are not in such good agreement, 
but differ by less than 2 standard deviations from the observed values. The value of 
Mb* — Mbo is 28% lower than experimental value, but it agrees with the result from a 
large and fine lattice[15]. 

For c-fiavored hadrons we take the lattice spacing determined by the charmonium 
system: 

= 0.75GeV (4.10) 

and obtain 

Md* - Mdo = 67(5)MeV, Md - Moo = 103(9)MeV, 

(4.11) 

Md* - Md, = 66(14)MeV, Ma, - M^o = 471(6)MeV. 
The experimental results are: 

Md* - Mdo = 143MeV, Md - Mdo = 104MeV, 

(4.12) 

Md* - Md, = 144MeV, Ma, - M^o = 420MeV. 

Our value for Md, — Mdo agrees well with the experimental value, but the predicted 
Ma, — Mdo is 12% larger the experimental value. More worse are the spin-splittings 
which are much lower than experimental values. Such discrepancies may be expected as 
we used lattice NRQCD upto order of 1/Mq. For c-quarks the effect from higher orders is 
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very significant because the charm mass is not very large. A recent study of quarkonium 
systems[16] also shows that the effect of higher order terms in 1/Mq is large. We also 
have calculated masses of c-flavored hadrons at /? = 6.8 with similar results. It shows that 
the action for lattice NRQCD employed here is not accurate enough for calculations of the 
spin-splitting for c-flavored hadrons. 

To obtain the absolute mass of a b- or c-flavored hadron we need to know Ag in 
Eq.(4.3). This quantity can be calculated perturbatively. It can also be extracted by 
studying the mass spectrum of the quarkonium system. We extracted this number from 
our previous study at /5 = 7.4: 

Ac = 1.59(11), A6 = 5.25(18). (4.13) 

It should be noted that effects from the spin-dependent interaction are neglected in the 
extraction. Adding the spin-dependent interaction has little effect on these numbers, and 
also these numbers have a quite large uncertainty. Using Eq.(4.13) we obtain 

Mdo = 1.79(8)GeV, Mgo = 5.20(15)GeV (4.14) 

where the lattice spacings in Eq.(4.10) and in Eq.(4.7) are used respectively. These pre- 
dictions are in good agreement with the experimental values: 

Mdo = 1.864GeV, M^o = 5.278GeV. (4.15) 

However, the accuracy of our predictions is not good because of the large statistical error 
of Aq. 

5. Summary 

In this work we used improved actions for glue and quarks to calculate masses of 
light, strange, c- and b- flavored hadrons, where for c- and b-quarks we employed lattice 
NRQCD. The actions for quarks are improved at tree-level to remove the effect at order 
of O(a^), while the gluonic action is improved at one-loop level. Tadpole improvement 
is implemented. The results obtained in this work for light and strange mesons are in 
agreement with experimental results. We obtain the mass spectrum of b-flavored hadrons 
with improved actions on a coarse lattice with the lattice spacing as 0.23fm, and it is in 
agreement with experiment. For c-flavored hadrons the large effect from higher orders 
neglected in lattice NRQCD prevents us from obtaining spin-splitting in mass comparable 
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with experiment. In this work we have shown that improved actions with tadpole improve- 
ment can be used not only in light hadron sectors and in quarkonium systems as shown 
already in previous studies, but also works in the case of hadrons containing one heavy 
quark. 

Acknowledgment: 

This work is supported by Australia Research Council. We thank our computer man- 
ager Dr. M. Munro for help with the implementation of our calculations. 



12 



Reference 



[I] G.P. Lepage and P.B. Mackenzie, Phys. Rev. D48 (1993) 2250 

[2] K. Symanzik in Mathematical Problems in Theoretical Phyics, ed. R. Schrader et al. 

Lecture Notes in Physics 153 (Springer Berlin 1982) 

Nucl. Phys. B226 (1983) 187 
[3] M. Alford et. al. Phys. Lett. B361 (1995) 87 
[4] J.P. Ma and B.H.J. McKellar, Melbourne Preprint UM-P-96/55 
[5] M. Alford, T. Klassen and G.P. Lepage, Nucl. Phys.(Proc. Suppl.) B47 (1996) 370 
[6] H.R. Fiebig and R.M. Woloshyn, Phys. Lett. B385 (1996) 273 
[7] H. Hamer and CM. Wu, Phys. Lett. B133 (1983) 351 

W. Wetzel, Phys. Lett. B136 (1984) 407 

T. Eguchi and N. Kawamoto, Nucl. Phys. B237 (1984) 609 
[8] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259 (1985) 572 
[9] M. Liischer and P. Weisz, Comm. Math. Phys. 97 (1985) 59 
[10] N. Cabibbo and E. Marinari, Phys. Lett. B119 (1982) 387 

[II] A. Prommer et.al. Int. J. Mod. Phys. C5 (1994) 1073 
[12] G.P. Lepage et al, Phys. Rev. D46 (1992) 4052 

[13] C.T.H. Davies et al, Phys. Rev. D50 (1994) 6963 

[14] P. Lacock and C. Michael, Phys. Rev. D52 (1995) 5213 

[15] A. Ah Khan et al., Phys. Rev. D53 (1996) 6433 

[16] H.D. Trottier, Preprint SFU HEP-131-96, |hep-lat 9611026 



13 



Figure Caption 



Fig.lA: The propagator C7r(t) with m = —.65. The points with error are measured, the 

hne is from the fit. The x-axis is for t. 

Fig. IB: The propagator Cp{t) with m = —.65. 

Fig.lC: The effective mass plot for the proton propagator with m = —.65, the effective 
mass is defined as mefi = ln(Cp(t — 1)/C{t)). The x-axis is for t. 

Fig. 2: The chiral extrapolation for p and proton. The upper line with data points is for 

proton. The lower line with data points is for p. The x-axis is read to be m^. 

Fig. 3: The chiral extrapolation for strange meson. 

Fig.4A: The effective mass plot for B with m = —0.65. 

Fig.4B: The effective mass plot for B* with m = —0.65. 

Fig.4C: The effective mass plot for A;, with m = —0.65. 

Fig.5A: The chiral extrapolation for b-flavored hadrons. The upper line with data points 
is for Aft, the middle one with data points is for B*, the lower one with data points is for 
B. The X-axis is read as m. 

Fig.5B: The chiral extrapolation for c-flavored hadrons. The upper line with data points 
is for Ac, the middle one with data points is for D*, the lower one with data points is for 
D. The X-axis is read as m. 
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